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ABSTRACT 


A  tactical  supercritical  flow  research  model  was  designed 
to  achieve  high  aerodynamic  efficiency  at  Mach  number  0.90  and  lift  coefficient 
0.50  using  numerical  transonic  design  procedures.  A  wind  tunnel  model  was 
built  and  tasted  In  the  Rockwell  International  Trl sonic  wind  tunnel  at  Mach 
numbers  of  0.40,  0.85,  0.90,  and  0.95.  Lift-drag  ratio  at  the  design  condition 
compared  favorably  with  theoretical  upper  bound  levels  for  the  test  arrangement 
and  conditions.  Six  component  force  measurements  for  the  joined  wing  are  well 
behaved  and  near  linear  over  an  angle-of- attack  range  of  -4<  a  <  10  degrees. 
Stable  stall  Is  naturally  achieved  for  the  arrangement.  Weak  changes  in  longi¬ 
tudinal  and  lateral -directional  stability  and  control  with  pitch  angle  exist. 
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NOMENCLATURE 


General  Symbols 
a.c. 

AR 

b 

8.P. 

c 

?,  MAC 

CG 

CP 

C“o 

So  »r  Sx 

Vl2 


X  =  0 

CMC 
S  =  £ 


V 

dCM/dCL 

e100 

FP 

F.S. 


Aerodynamic  center 
2 

Aspect  ratio,  b  /Srgf 
Wing  span 

Buttock  plane,  inches 
Surface  chord 
Mean  aerodynamic  chord 
Center  of  gravity 
Center  of  pressure 

0  Percent  suction  drag  or  drag  at  CL  ~  0 
100  Percent  suction  or  vortex  drag 
Drag  due  to  lift  factor 
Sectional  lift  coefficient,  l/qc 
Lift  curve  slope,  dC^/do,  per  degree 
Lift  coefficient  at  a  =  0  degrees 
Pitching  moment  slope,  dC^do,  per  degree 
Pitching  moment  at  a  =  0  degrees 
Pitching  moment  at  C^  =  0 
Sectional  normal  force  coefficient 
Pressure  coefficient,  (P-P*^i 
Critical  pressure  coefficient 
Static  longitudinal  stability 
Trefftz  plane  spanload  efficiency 
Full  potential 
Fuselage  station,  inches 
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RN,  R. 


X.Y.Z 


Wing  Incidence,  negative  nose  down 
Lift  to  drag  ratio,  C^/Cq, 

Mach  number 

Local  upper  surface  Mach  number 

Static  pressure 

Dynamic  pressure,  psf 

Total  velocity 

Reynolds  number  per  foot 

Reynolds  number  based  on  mean  aerodynamic  chord 
2 

Reference  area,  ft 

Aft  wing  reference  projected  trapezoidal  wing  area 

Forward  wing  reference  projected  trapezoidal  area 

Wing  thickness  ratio 

Transonic  small  disturbance 

Water  plane,  Inches 

Wing  station,  Inches 

Cartesian  body  axis  axial,  lateral,  vertical 
coordl nates 


Angle-of-attack,  degrees 
Angle-of-sldeslip,  degrees 
Increment 

Surface  Deflection,  perpendicular  to  hinge 
line,  degrees 

Oownwash  angle,  degrees 

Fraction  of  span,  2y/b 

Sweep  angle,  degrees 
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NOMENCLATURE  (CONTINUED) 


6 


Twist  angle,  positive  nose  up 


¥ 


Taper  ratio 

Yaw  angle,  (-0),  degrees 


Coefficients 

Cg  Chord  force  coefficient 

Cq  Drag  Coefficient, 


Chord  Force 
qS 


CL 

C£ 

CM 

CN 

Cn 

CP 

cy 


Lift  coefficient, 

Rolling  moment  coefficient. 

Pitching  moment  coefficient, 

Normal  force  coefficient, 

Yawing  moment  coefficient, 
Pressure  coefficient, 

Side  force  coefficient, 


Lift 

IT 

Rolling  Moment 

qSb 

Pitching  Moment 


mg 

j5c 


Normal  Force 
- q5 - 

Yawing  Moment 
q5E 

P-P„ 


Side  Force 
- qS - 


Subscripts 

A 

AL 

B 


Aft  wing  panel 
Attachment  line 
Body  axes 
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INTRODUCTION 
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The  joined  wing  is  a  general  concept  which  combines  two  wings,  a  fuselage 
and  a  vertical  fin.  The  wings  form  a  diamond  shape  in  both  planform  and  front 
view.  A  limited  amount  of  low  speed  design,  wind  tunnel  test  and  evaluation, 
and  several  structural  studies  are  discussed  in  reference  1  and  in  references 
cited  therein.  Joined  wing  aerodynamic  technology  offers  potential  for 
improved  induced  drag  efficiency  for  transonic  flight.  An  aerodynamic  study,  of 
the  application  of  the  joined  wing  concept  to  a  fighter  can  be  used  to  explore 
this  potential  and  to  define  the  unique  characteristics  of  a  joined  wing 
arrangement.  A  Mach  number  of  0.90  lift  coefficient  of  0.50  design  point 
was  selected  to  focus  the  study  and  corresponds  to  a  condition  between  cruise 
and  high  g  maneuvers.  Contour  definitions  by  numerical  transonic  design,  and 
experimental  validation  by  subsonic/transonic  wind  tunnel  test,  are  required  to 
validate  the  highly-nonplanar,  strongly-interacting  generic  joined  wing 
arrangement. 

Configuration  aerodynamic  development,  numerical  wing  design  and  experi¬ 
mental  test  validation  are  described  in  the  text  of  the  report.  Specific 
methodology  used,  viscous  drag  analysis  and  model  design  and  construction 
details  are  presented  in  appendices  A,  8,  and  C,  respectively. 
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Section  II 

CONFIGURATION  DEFINITION 


•Representative  geometric  data  from  current  fighters  (e.g. ,  body  length 
to  wing  span,  body  width  to  wing  span,  tall  volume  etc.)  were  used  to  define 
a  reference  monoplane  and  a  jolned-wing  configuration.  A  common  body  was  used 
for  both  configurations.  Standard  longitudinal  locations  of  the  aerodynamic 
surfaces  were  used  for  the  monoplane. 

The  joined  wing  configuration  was  selected  to  have  the  same  span  as  the 
reference  monoplane  and  a  projected  total  (i.e. ,  to  the  plane  of  symmetry) 
plan  aria  of  both  panels  equal  to  the  total  wing  area  plus  exposed  horizontal 
tall  area  of  the  reference  monoplane.  The  joined  wing  panel  leading  edge 
sweep  was  nominally  selected  to  be  ±45  deg;  ,  based  on  transonic  perpen¬ 
dicular  sectional  lift  coefficients  (C£  -  *  cos2Ac/t2)  and  supersonic  wave 

drag  considerations.  Ten  degrees  of  dihedral  and  20  degrees  of  anhedral  were 
used  for  the  forward  and  aft  wing  panels,  respectively,  to  achieve  substantial 
reductions  of  vortex  drag  due  to  nonplanar  effects.  The  vertical  tall  size 
was  a  geometric  fallout  of  the  30  degrees  Included  angle  of  the  joined  wing. 
Nominal  thickness  ratio  t/c  a  0.05  was  selected  based  on  current  fighter  design 
practice.  Configuration  arrangement  of  the  reference  monoplane  and  the  joined 
wing  are  shown  on  figures  1  and  2,  respectively.  Note  that  the  dimensions 
shown  correspond  to  a  l/10th  scale.  The  wind  tunnel  model  tested  was  a  l/15th 
scale  as  described  later. 

Arrangement  considerations  for  the  joined  wing  Included  a  level  of  low 
speed  static  longitudinal  Instability  and  lateral -directional  stability  that 
were  similar  to  the  reference  monoplane  and  consistent  with  current  tactical 
advanced  design  studies.  A  comparison  of  the  various  linear  aerodynamic 
characteristics  computed  by  the  method  of  reference  2  are  presented  on  table  I. 

The  pitch  damping  derivative,  ,  is  typically  twice  that  of  the  reference 
monoplane.  q 
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Section  III 
JOINED  WING  DESIGN 


The  numerical  aerodynamic  design  methodology  used  for  the  present  study 
is  summarized  In  appendix  A. 


LINEAR  INITIALIZATION 

Two  linear  nonplanar  finite  element  models  were  used'  for  the  study. 
The  analysis  model  consists  of  a  slender  body  (not  shown).  Interference  shell, 
forward  and  rear  wings,  and  a  vertical  tall  as  shown  In  figure  3.  Vortex  drag 
minimization  (reference  3)  uses  a  flat  plate  fuselage  and  forward  and  aft-wings 
as  shown  In  figure  4.  The  vertical  tail  Is  not  simulated  because  It  has  no 
Influence  on  the  optimization.  In  order  to  maintain  solution  consistency 
between  models,  10  uniform  chordwlse  by  10  uniform  spanwlse  panels  were  used 
for  the  forward  and  aft  wings. 

Linear  vortex  drag  optimizations  were  performed  to  satisfy  necessary 
(far  field)  design  objectives  at  M  =  0.90  and  C,  *  0.50.  The  following 
constraints  were  considered  to  Insure  a  practical  resftilt: 

1.  CM  =  0  at  the  design  CL 

2.  High  lift  stable  pitching  moment  break 

3.  Moderate  supersonic  longitudinal  stability  levels 

4.  Moderate  twist 

5.  Sectional  loadings  that  are  consistent  with  moderate  upper 

surface  supercritical  Mach  numbers  (M^  <  1.3) 

Equal  joined  wing  projected  panel  areas  were  considered  initially  (cases  1 
through  4  of  table  II).  The  Impact  of  longitudinal  balance  on  trimmed  span 
load  efficiency  and  the  relative  peak  sectional  lift  (cases  1  through  3)  are 
presented  on  figures  5  and  6  respectively.  As  shown,  unstable  balance  Is 
Indicated  which  was  compatible  with  low  speed  stable  stall  (since  the  forward 
wing  panel  progressively  separates  from  tip  to  root  and  Is  located  ahead  of 
the  center  of  gravity)  and  moderate  supersonic  longitudinal  stability. 
Attempts  to  more  uniformly  balance  the  loading  between  the  forward  and  aft 
panels  (case  4),  In  the  Interest  of  reducing  the  supersonic  upper  surface  Mach 
number  on  the  forward  wing,  resulted  in  unacceptable  span  load  efficiency 
penalties  (figure  5)  and  excessive  twist.  Two  alternative  approaches  were 
examined  to  reduce  the  loading  on  the  forward  panel: 

1.  Reduce  the  sectional  design  lift  coefficient  (which  is  uniquely 
related  to  the  gross  lift  coefficient  for  a  given  optimum 
solution). 

2.  Redistribute  wing  area  to  60  percent  on  the  forward  panel  and 
40  percent  to  the  aft  panel  (case  5). 
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The  results  of  the  latter  trade  are  presented  on  figure  5  and  indicate 
penalties  on  span  load  efficiency  relative  to  the  equal  area  case.  The 
reduced  design  C,  option  was  consequently  selected.  A  local  upper  surface 
Mach  number  of  l1.?5  on  the  forward  panel  was  tentatively  selected  and 
corresponds  to  a  gross  design  C,  =  0.45  for  the  20  percent  unstable  M  =  0.90 
balance  (case  6).  Table  II  presents  a  summary  of  the  initial  optimization 
cases  considered. 

The  design  C,  was  subsequently  increased  to  0.50  to  reflect  a  constant 
more  positive  lower  surface  candidate  pressure  than  originally  assumed  (see 
figure  7).  As  indicated,  this  revision  allows  an  increase  in  sectional  lift 
coefficient  for  a  fixed  upper  surface  Mach  number. 

A  series  of  linear  vortex  drag  optimizations  were  subsequently  performed 
at  M  -  0.90,  C,  =  0.50  with  a  longitudinal  instability  of  Q.2QF.  Several 
design  refinements  were  incorporated  subject  to  maintaining  the  previously 
derived  potential  drag  due-to-lift.  They  were:  smooth  variation  of  spanwise 
twist,  low  toot  chord  incidence,  and  supercritical  type  camber  on  the  forward 
surface.  The  resulting  linear  design  is  presented  on  figure  0.  The  impact  of 
the  previously  cited  constraints  on  span  load  efficiency  was  negligible. 


LINEAR  INDUCTION  EFFECTS  AND  EFFECTIVE  TWIST 


The  method  used  to  account  for  the  wing-wing  interaction  in  the  Transonic 
Small  Disturbance  (TSD)  analysis  is  an  immersion  or  apparent  twist  philosophy. 
The  upwash/downwash  distribution  of  the  aft  wing  on  the  front  and  the  front 
wing  on  tie  aft  are  calculated  with  linear  theory  and  then  applied  as  twist 
Increments  to  the  geometry  used  in  the  TSD  code. 

The  induction  effects  are  calculated  for  the  candidate  linear  design  case 
using  a  mixed  boundary  condition  analysis  (reference  3).  The  process  Imposes  a 
r 

Pnet  *  0  constra'*nt  on  the  surface  where  the  downwash  effects  are  to  be 

r 

defined,  and  maintains  the  design  Pne^is  on  the  other  surface(s).  The  induc¬ 
tion  effects  of  the  aft  wing  on  the  front  wing  are  negligible.  The  effects  of 
the  front  on  the  rear  wing  are  given  on  figure  9.  Induced  camber  is  negligible 
for  either  panel. 


DESIGN  PRESSURE  DISTRIBUTION 


Evaluation  of  M  =  0.90  candidate  sectional  supercritical  pressure  distri¬ 
bution  for  the  forward  and  aft  wing  panels  was  performed  for  the  nominal  model 
test  Reynolds  number.  The  design  sectional  lift  coefficients  were  0.60  and 
0.33  for  the  forward  and  aft  panels  respectively,  for  a  design  gross  lift 
coefficient  of  0.50.  These  pressure  distributions  correspond  to  midspan  values 
resulting  from  minimization  of  the  trimmed  vortex  drag  for  the  full  three 
dimensional  arrangement.  Flat  top  upper  surface  distributions  with  a  linear 
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recompression  beginning  at  60  percent  chord  were  initially  considered.  Local 
Mach  numbers  of  1.19  for  the  swept  aft  wing  and  1.07  for  the  swept  forward 
panel  resulted  for  the  sectional  C.  of  Interest  considering  representative 
lower  surface  pressure  distribution* levels.  The  resulting  candidate  design 
pressure  distributions  are  presented  on  figures  10  and  11.  Local  sonic  and 
leading  edge  attachment  line  (1.e.,QAL  =  Q*  sinALE)  conditions  are  indicated. 

An  integral  laminar/turbulent  yawed  wing  boundary  layer  analysis 
(reference  4)  was  performed  for  a  joined  wing  mean  aerodynamic  chord  Reynolds 
number  of  4.1  x  10®  which  corresponds  to  a  minimum  dynamic  pressure  tunnel 
operating  condition  to  reduce  model  stress  levels.  Transition  was  enforced  at 
approximately  3.5  percent  chord.  The  results  of  figure  12  indicate  the  flow 
remains  fully  attached  for  the  forward  surface  with  separation  occurring  down¬ 
stream  of  95  percent  chord  on  the  upper  surface  of  the  aft  wing.  Although  the 
latter  panel  is  more  lightly  loaded,  its  higher  trailing  edge  sweep  offsets  the 
reduction  in  adverse  pressure  gradient. 

The  results  of  the  viscous  analysis  indicate  the  candidate  design  pressure 
distributions  of  figures  10  and  11  satisfy  the  sufficient  condition  of  good 
flow  quality  considering  real  fluid  effects.  The  potential  far  field  optimi¬ 
zation  results  of  figure  3  satisfy  the  necessary  design  conditions.  Taken 
together,  both  necessary  and  sufficient  considerations  are  met  and  provide 
satisfactory  initialization  of  the  three  dimensional  transonic  design.  The 
boundary  layer  displacement  thickness  growth  for  the  upper  and  lower  surface  of 
figure  12  indicate  a  candidate  airfoil  should  have  a  nominal  trailing  edge 
thickness  of  1.5  percent  chord  to  permit  undercutting  of  the  design  contour  to 
minimize  weak  viscous  Interaction  effects  on  the  transonic  potential  design. 


CANDIDATE  SUPERCRITICAL  AIRFOILS 

The  nominal  airfoil  design  conditions  for  M  -  0.90  C.  =  0.50  trimmed 
condition  with  a  20  percent  longitudinal  unstable  balance  and  a  thickness 
ratio  of  5  percent  are  summarized  on  table  III.  The  candidate  streamwise  air¬ 
foil  sections  for  the  forward  and  aft  swept  wings  are  presented  in  figures  13 
and  14  respectively.  These  contours  in  conjunction  with  the  effective  twist 
distributions  of  figures  8  and  9  are  used  to  initialize  the  three  dimensional 
potential  transonic  small  disturbance  design  of  the  aft  swept  wing-body  and 
forward  swept  wing  vertical  tail  fairing  described  in  the  next  section. 


TRANSONIC  SMALL  DISTURBANCE  DESIGN 
FORWARD  WING 

The  transonic  small  disturbance  (TSD)  model  for  the  front  wing-body 
combination  Is  shown  on  figure  15.  Figure  15a  shows  the  wing  planform  and  the 
prism  used  to  satisfy  the  body  boundary  conditions.  Figure  15b  shows  the 
fuselage  modeling.  The  forward  wing  was  initially  modeled  with  the  candidate 
airfoil  of  figure  13  across  the  entire  span  in  conjunction  with  the  lesser 
optimum  twist  of  figure  8.  The  resulting  analysis  (reference  5)  is  shown  on 
figure  16. 
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To  reduce  the  aft  loading  of  the  Inboard  airfoils  the  camber  was  progressively 
decreased.  The  Inboard  span  stations  have  a  relatively  strong  shock  at  90 
percent  chord  which  required  further  design  work  to  weaken.  The  mid-span 
region  developed  the  desired  flat  top  pressure  distribution  although  terminal 
upper  surface  adverse  pressure  gradients  exceed  the  candidate  distribution  and 
required  moderation.  The  outboard  region  meets  the  design  goal  and  required 
no  additional  refinements. 

Further  design  effort  attempted  to  eliminate  the  Inboard  trailing  edge 
shock.  The  target  upper  surface  pressure  coefficient  was  specified  aft  of 
20  percent  chord  at  each  of  the  Inboard  TSD  defining  stations.  The  resulting 
airfoils  were  unrealistic.  Subsequent  effort  was  directed  at  more  modest  root 
shock  weakening  movement  of  the  shock  towards  the  trailing  edge.  A  candidate 
5  percent  airfoil  was  derived  from  Inverse  analysis  (reference  6)  In  conjunc¬ 
tion  with  twist  adjustments  to  approximate  the  linear  optimum  sectional  normal 
force  distribution.  Comparison  of  the  results  indicate: 

1.  The  strength  of  the  Inboard  shock  and  adverse  pressure 

gradient  has  been  moderately  reduced. 

2.  The  Inboard  change  of  airfoil  caused  a  degradation  of  the 

flow  quality  In  the  midspan  region  of  the  wing. 

3.  The  outboard  section  of  the  wing  was  unchanged. 

The  Impact  on  the  center  section  was  not  considered  serious  since  subsequent 

design  effort  was  to  be  directed  at  reducing  the  trailing  edge  recompression 

problems  In  this  region. 

The  midspan  refinement  effort  concentrated  on  the  area  from  47  to  61  percent 
of  the  span.  An  Inverse  design  was  performed  for  the  target  upper  surface 
pressure  distribution  defined  from  approximately  1  percent  chord  to  the  trail¬ 
ing  edge.  This  solution  defined  an  airfoil  that  exhibited  Improved  center 
section  flow  quality.  At  this  point  In  the  design,  the  model  stress  analysis 
Indicated  the  wing  tip  area  required  an  Increase  In  thickness  from  5  to  7 
percent  with  a  linear  variation  to  5  percent  at  the  midspan.  Subsequent  analysis 
of  this  modification  Indicated  only  a  minor  Impact. 

A  three  dimensional  finite  difference  boundary  layer  analysis  (reference  7) 
was  performed  at  the  tunnel  test  mean  aerodynamic  chord  Reynolds  number  of 
4.1  x  10°  to  verify  fully  attached  flow  would  be  achieved  except  inboard 
near  the  trailing  edge  on  the  upper  surface.  The  airfoils  resulting  from  the 
potential  design  were  undercut  by  this  boundary  layer  displacement  thickness 
to  minimize  weak  viscous  Interaction  effects.  The  resulting  trailing  edge 
In  the  outboard  region  had  insufficient  structural  thickness.  The  lower  sur¬ 
face  was  subsequently  modified  for  x/c  >  80  percent.  The  Increased  trailing 
edge  thickness  had  the  additional  benefit  of  reducing  the  adverse  pressure 
gradients  on  the  lower  surface.  The  final  forward  wing  transonic  design 
geometry  Is  shown  on  figure  17.  The  design  flow  quality  is  presented  on 
figure  18.  Figures  19  and  20  show  the  boundary  layer  displacement  thickness 
and  undercut  airfoils,  respectively. 
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The  aft  wing  was  initially  modeled  as  a  trapezoidal  planform  with  the 
candidate  airfoil  of  figure  14  across  the  entire  span.  The  initial  twist  is 
given  in  figure  9.  As  a  result  of  not  achieving  the  desired  flow  quality  when 
matching  the  linear  theory  optimum  span  load  shape,  several  modifications  were 
made  to  the  camber,  twist  and  planform.  The  camber  was  decreased  at  the  root 
and  the  tip  to  increase  the  forward  loading  of  the  lower  surface.  The  twist 
was  also  adjusted  to  give  the  desired  sectional  lift  distribution  across  the 
span.  The  revised  planform  incorporates  an  inboard  glove  to  reduce  the 
strength  of  the  flow  singularity  at  the  root,  and  a  central  body  fairing  of 
circular  cross  section  was  added  due  to  model  structural  wing- tail  attachment 
considerations,  figures  21a  and  21b.  The  status  airfoil  at  n  =  0.04  (the  side 
of  the  fairing)  is  presented  in  figure  21c.  The  section  varies  linearly  to  the 
airfoil  at  n  *  0.25.  The  airfoil  of  figure  21d  is  used  from  n  a  0.25  to 
0.80.  The  airfoil  then  varies  linearly  to  that  of  figure  21e  at  the  tip. 
Status  geometric  twist  is  compared  with  the  linear  optimum  on  figure  21f. 
Transonic  small  disturbance  (reference  5)  span  loading  and  upper  and  lower 
surface  pressure  distribution  Isometrics  are  given  on  figures  21g  and  21h 
respectively.  Figures  211  through  21k  show  sectional  pressure  distributions 
at  three  representative  span  stations.  The  inboard  span  stations  still 
Indicated  a  shock  while  the  midspan  region  develops  the  desired  flat  top 
pressure  distribution  with  a  mild  leading  edge  peak. 

Further  design  effort  was  directed  at  eliminating  the  root  shock, 
reducing  the  leading  edge  peak  pressure  in  the  center  section,  and  increasing 
the  upper  surface  leading  edge  pressure  at  the  tip.  A  transonic  small  distur¬ 
bance  Inverse  (reference  6)  analysis  was  employed  to  determine  changes  to  the 
airfoil  geometry  to  eliminate  the  root  shock  and  used  to  reduce  the  upper 
surface  pressure  peak  at* the  leading  edge  of  the  wing  center  section.  Tip 
airfoil  redesign  was  directed  at  increasing  the  loading  In  the  leading  edge 
region.  Since  further  Increase  in  twist  was  undesirable  based  on  tip  fairing/ 
attachment  considerations,  effort  concentrated  on  decambering  the  airfoil 
section  while  maintaining  the  linear  optimum  sectional  normal  force 
coefficient  distribution. 

A  three  dimensional  finite  difference  boundary  layer  analysis 
(reference  7)  was  gthen  performed  at  a  test  mean  aerodynamic  chord  Reynolds 
number  of  4.1  x  10°  to  verify  fully  attached  flow  was  essentially  achieved  and 
define  associated  boundary  layer  displacement  thickness  distribution  for  the 
upper  and  lower  surfaces.  Separation  occurs  nominally  in  the  last  5  percent  on 
the  upper  surface  and  is  fully  attached  on  the  lower  surface.  The  inviscid 
design  airfoils  are  corrected  for  weak  viscous  interaction  effects  by  under¬ 
cutting  the  coordinates  by  the  displacement  thickness. 

A  7  percent  thick  tip  section  was  subsequently  investigated  due  to  model 
pressure  instrumentation/structural  considerations.  The  design  flow  quality 
was  not  significantly  affected.  The  final  design  geometry  incorporating  a 
7  percent  thick  tip  section  linearly  decreasing  to  a  base  5  percent  airfoil 
at  50  percent  semispan  is  summarized  in  figure  22  and  23.  Figure  24  presents 
the  final  design  loading  and  flow  quality.  Figures  25  and  26  present  the  aft 
wing  boundary  layer  displacement  thickness  and  undercut  airfoils 
respectively. 
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THEORETICAL  DESIGN  PERFORMANCE 

A  theoretical  figure  of  merit  for  the  joined  wing  will  now  be  defined. 
Specifically,  it  is  the  highest  (or  ultimate)  L/D  that  would  be  expected  at  the 
design  condition.  The  evaluation  includes  fully  turbulent  skin  friction, 

profile  and  potential  wake  losses  but  does  not  account  for  viscous  component 
Interference  or  separation  drag.  As  such  the  performances  is  upper  bound  in 
nature.  Skin  friction  drag  is  evaluated  (reference  2)  for  a  flat  plate  condi¬ 
tion  in  conjunction  with  a  component  build  up  approach.  The  result  (figure  27) 
is  corrected  to  account  for  the  design  pressure  distribution  and  sweep  using 
the  profile  drag  method  described  in  appendix  B.  Induced  drag  is  calculated 
using  trimmed  potential  theory  lifting  efficiency.  The  drag  is  then  the  sum 

c  c  c 

of  skin  friction  (  Df),  profile  drag  (  Dp)  and  the  induced  (  Dj)  drag. 

The  constituent  drag  pieces  and  the  trimmed  upper  bound  L/D  at  the  design 
lift  coefficient  are  summarized  in  figure  28. 


TRANSONIC  REDESIGN 


The  results  of  the  first  wind  tunnel  test  are  presented  in  reference  8, 
appendices  C,  D,  and  E.  Pertinent  design  results  are  shown  in  figure  29.  Lift 
curve  slope  and  longitudinal  stability  were  as  expected  while  the  intercepts 
Cl  ,|.J  Cy 

Lo  ana  no  were  not  matched  by  linear  analysis  for  the  design  twist  and 
as  built  camber.  As  a  consequence  self-trim  was  not  achieved. 

Examination  of  the  wind  tunnel  pressure  data  revealed  the  forward  wing  was 
more  heavily  loaded  and  the  aft  wing  more  lightly  loaded  than  Intended.  This 
resulted  In  the  forward  wing  developing  a  swept  shock  (figure  30)  at  the  design 
C,.  The  front  wing  center  section  was  intended  to  operate  at  lower  angle-of- 
attack  and  the  rear  wing  airfoils  at  higher  angle-of-attack  due  to  the  elastic 
twist  increments  (figure  31).  When  structurally  bench  tested,  (see  appendix 
C)  the  model  was  not  nearly  as  flexible  as  a  beam  structural  analysis  had 
predicted. 

Calculations  were  subsequently  performed  for  the  as  built  (i.e.  ,  jig) 
twist.  The  linear  spanloading  (figure  32)  showed  this  was  indeed  the  cause 
of  the  loading  anomalies  observed.  The  jig  twist  was  then  analyzed  using 
transonic  small  distribution  theory  at  the  wind  tunnel  angle-of-attack  for 
the  design  lift  coefficient.  This  result  (figure  33)  produced  the  anomalous 
results  that  no  shocks  would  be  anticipated  for  either  the  design  or  jig  twist 
on  the  forward  wing.  Conservative  full  potential  analysis  (reference  9)  was 
subsequently  performed. 

These  results  (figure  34)  Indicate  an  essentially  shockless  condition  on 
the  forward  wing  with  the  design  twist  and  angle-of-attack  and  a  shocked  flow 
for  the  jig  twist  and  test  angle-of-attack  for  C,  =  0.50.  A  comparison  of  the 
forward  wing  numerical  shock  capturing  differences  between  transonic  small 
disturbances  and  full  potential  theory  at  the  design  C,  =  0.50  for  the  as 
built  twist  are  summarized  on  figure  35.  The  latter  analysis  was  consequently 
used  for  redesign  effort.  The  aft  wing  flow  (figures  36  and  37)  was  adequately 
predicted  by  TSD  and  verified  to  be  shockless  at  the  design  condition. 
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Differences  between  design  and  jig  twist  for  the  aft  wing  are  offset  by 
the  Increased  angle-of-attack  required  to  achieve  C,  =  0.50.  The  opposite  is 
true  for  the  forward  wing.  This  consideration  in  conjunction  with  substantially 
higher  design  loading  accounts  for  its  elastic  sensitivity. 

Since  the  twist  of  the  high  strength  steel  wings  could  not  be  readily 
changed,  the  use  of  negative  wing  Incidence  to  reduce  the  forward  panel  load¬ 
ing  was  examined  as  an  alternative.  Linear  theory  indicated  a  1  degree  nose 
down  Incidence  on  the  forward  wing  and  a  0.25  degree  nose  up  change  on  the  aft 
wing  would  approximate  the  design  span  load  shape.  The  resources  available 
did  not  permit  Incorporation  of  the  aft  wing  modification.  The  linear  span  load 
for  a  -I  degree  change  on  the  forward  wing  Is  presented  on' figure  38.  The 
predicted  longitudinal  results  are  presented  on  figure  39.  The  self  trim  condi¬ 
tion  at  the  design  C.  Is  essentially  realized.  Full  potential  analysis  of  the 
modified  design  Is  presented  on  figure  40  and  Indicates  an  essentially 
shockless  flow. 

The  results  of  the  second  wind  tunnel  test  are  presented  In  appendices  G, 
H,  and  I  of  reference  8.  Pertinent  design  results  are  shown  In  figure  41. 
Lift  curve  slope  and  longitudinal  stability  were  unchanged  from  previous 
results.  The  Intercepts,  L  and  uMfl,  again  were  not  matched  by  linear 
analysis  for  the  jig  twist,  wing  Incidence  change,  and  as  built  cambers.  The 
design  Is  considerably  closer  to  achieving  self  trim.  As  before,  a  swept  shock 
exists  on  the  forward  wing  but  appears  to  be  weaker  for  the  center  section  due 
to  the  forward  wing  still  being  more  heavily  loaded  than  desired.  A  detailed 
discussion  of  the  second  test  entry  surface  pressure  measurements  and 
associated  oil  flow  data  Is  presented  In  section  VI.  Further  flow  quality 
progress  (i.e. ,  shock  weakening)  requires  achieving  the  design  twist  on  both 
the  forward  and  aft  panels  to  realize  the  Intended  loading  distribution  between 
the  two  Surfaces. 

The  aft  wing  flow  quality  Is  good  and  consequently  does  not  require 
further  attention.  This  Is  due  In  part  to  the  panel  being  more  lightly 
loaded  for  the  design  condition. 
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Section  IV 

REFERENCE  MONOPLANE  DESIGN 


The  reference  monoplane  of  section  II  was  derived  from  representative 
geometric  data  of  current  fighter  aircraft.  The  resulting  arrangement  Is  shown 
on  figure  1.  Representative  geometric  properties  Included:  body  length  and 
width  to  wing  span,  and  horizontal  and  vertical  tail  volumes. 


LINEAR  INITIALIZATION 

Necessary  monoplane  lifting  conditions  were  satisfied  by  optimizing  the 
far  field  vortex  drag  at  H  *  0.90  and  C,  =  0.50.  The  drag  minimization  model 
uses  a  flat  plate  fuselage,  wing  and  horizontal  tall  components.  A  vertical 
tail  was  not  simulated  since  it  has  no  Influence.  The  arrangement  was  paneled 
In  a  similar  manner  as  the  joined  wing  with  10  uniform  chordwlse  and  spanwlse 
finite  elements  for  the  wing  (see  figure  42). 

The  following  constraints  were  Imposed  to  Insure  a  practical  result: 

1.  CM  =  0  at  the  design 

2.  Section  loadings  consistent  with  moderate  upper  surface 
supercritical  Mach  numbers  (H£  <  1.30) 

3.  No  fuselage  twist  or  camber 

4.  Horizontal  tall  deflection  only  (that  is,  no  twist  or  camber) 

The  effect  of  longitudinal  balance,  on  constrained  span  load  efficiency 
Is  presented  on  figure  43.  A  20  percent  unstable  balance  was  subsequently 
selected.  The  section  loadings  for  this  case  are  presented  on  figure  44. 


DESIGN  PRESSURE  DISTRIBUTION 


Evaluation  of  the  design  sectional  supercritical  pressure  distribution 
(figure  45)  was  performed  for  the  mean  chord  test  Reynolds  number  of  7.3  x  10°. 
The  reference  monoplane  operates  at  a  higher  Reynolds  number  due  to  a  longer 
chord  than  the  joined  wing.  The  sectional  lift  coefficient  for  the  wing  at  the 
design  point  Is  0.55.  This  pressure  distribution  corresponds  to  a  near  midspan 
value  resulting  from  the  minimization  of  trimmed  vortex  drag  and  is  comparable 
in  magnitude  and  shape  to  the  distribution  on  the  forward  panel  of  the  joined 
wing.  A  local  upper  surface  Mach  number  of  1.16  resulted. 
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A  laminar/turbulent  yawed  wing  boundary  layer  analysis  (reference  4)  was 
initially  performed  for  the  same  unit  Reynolds  number  and  a  wing  mean  sweep  of 
33.6  degrees  corresponding  to  the  same  leading  edge  sweep  as  the  joined  wing. 
A  revised  reference  monoplane  having  the  same  mean  sweep  was  subsequently 
analyzed  in  order  to  achieve  similar  transonic  perpendicular  sectional  lift 

coefficients  (C.  *  *  cos2  A^)  and  comparable  supersonic  wave  drag  which 

is  a  function  of  the  maximum  thickness  line  which  is  approximately  equal  to 
the'  mean  sweep  line.  Transition  was  enforced  at  approximately  3.S  percent 
chord  for  both  cases.  The  analysis  indicated  fully  attached  flow  would  be 
expected  for  the  design  pressure  distribution  of  figure  45  for  either 
planform.  Profile  drag  results  using  the  analysis  of  appendix  8  are  summarized 
below. 


Planform 

ale 

(deg) 

Ac/2 

(deg) 

Rer 

CDp 

CDp 

Initial 

45 

33.6 

7.3  x  106 

0.01098 

0. 01170 

Revised 

53.2 

45 

7.3  x  106 

0.01098 

0.01288 
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Section  V 

JOINED  WING/REFERENCE  MONOPLANE  COMPARISON 
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Theoretical  Induced  drag  and  fully  turbulent  profile  drag  were  calculated 
ait  M  ■  0.90  and  C,  ■  0.50  for  the  two  configurations  for  the  same  wind  tunnel 
test  unit  Reynolds  number.  The  results  are  summarized  In  figure  46a.  The 
joined  wing  has  an  11  percent  Induced  drag  advantage  compared  to  the  monoplane 
of  the  same  span  due  to  the  nonplanar  nature  of  the  arrangement.  Viscous 
comparison  Indicate  the  joined  wing  has  15  percent  higher  viscous  drag  than  the 
monoplane  as  a  result  of  12  percent  higher  surface  area  for  the  same  plan  area 
and  a  42  percent  lower  mean  chord  Reynolds  number.  Ultimate  L/D's  at  the  design 
C.  ■  0.50  for  the  joined  wing  and  reference  monoplane  are  12.97  and  13.24, 
respectively.  These  levels  neglect  viscous  component  Interference,  separation 
drag,  etc.,  and  consequently  are  upper  bound  In  nature. 

A  similar  calculation  was  performed  for  full  scale  conditions  as  a  function 
of  flight  altitude.  The  results  are  presented  on  figure  46b  and  Indicate  the 
joined  wing  and  reference  monoplane  have  essentially  equal  upper-bound  aerodynamic 
efficiency  between  sea  level  and  50,000  feet  altitude.  Fixed  altitude  compari¬ 
sons  correspond  to  constant  normal  load  factor  In  order  to  maintain  C,  ■  0.50.* 
Constant  Reynolds  number  comparison  corresponds  to  variable  load  factor  due  to 
the  change  In  altitude  required  to  operate  at  the  design  lift  coefficient. 

The  previous  results  do  not  weigh  more  extensive  passive  laminar  flow 
extents  for  the  joined  wing  relative  to  the  monoplane  that  may  result  from 
operating  at  lower  chord  Reynolds  number. 


*wing  area  variations  were  not  considered  In  the  construction  of  the  results  of 
figure  46b. 
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Section  VI 

JOINED  WING  COMPARISON  WITH  THEORY 
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FORCE  DATA 

The  chord  plane  finite  element  models  of  figures  3  and  4  were  used  to 
predict  the  longitudinal  and  lateral -directional  characteristics  of  the 
joined  wing. 


LONGITUDINAL 


Measured  lift  and  pitching  moment  are  slightly  nonlinear  in  the  angle-of- 
attack  range  -4  <  a<  +10  degrees.  The  stall  Is  stable  because  the  outboard  sec¬ 
tion  of  the  forward  wing  stalls  first  end  is  forward  of  the  moment  reference 
point.  These  characteristics,  and  the  effect  of  the  forward  wing  Incidence 
change  have  been  presented  earlier  In  figures  29  and  39.  The  measured 
longtudlnal  derivatives  of  lift  curve  slope,  ''L ,  and  pitching  moment  slope, 
dCg/dC,,  are  expected  to  agree  most  closely  with  results  of  linear  theory 
calculations  at  angle-of- attack  and  sideslip  near  a  -  0  and  p  =  0  degrees. 
These  measured  derivatives  and  lift  and  pitching  moment  at  zero  angle  of 


attack  are  compared  on  table  IV  with  linear  theory  calculations.  Measured 
lift  curve  slopes  exceed  linear  theory  calculation  by  3  to  9  percent.  Measured 
lift  coefficients  at  zero  angle  of  attack,  UL  ,  are  40  to  60  percent  lower 
than  predicted  by  linear  analysis  using  as-Built  cambers.  The  neglect  of 
boundary  layer  displacement  thickness  decambering  effects  for  supercritical 
type  pressure  distributions  Is  believed  to  be  responsible  for  this  discrep¬ 
ancy.  The  Intent  was  to  select  the  moment  reference  point  so  that  the  joined 
wing  model  Is  20  percent  unstable  at  M  =  0.90.  The  dCM/dC.  comparison  of 
table  IV  shows  that  the  moment  reference  point  should  mavff  been  selected 


4  percent  TT  further  forward.  Measured  pitching  moment  at  zero  angle-of-attack, 
UM I,  was  within  0.008  of  linear  theory  prediction.  The  boundary- layer  Induced 
lift  loss  is  In  the  same  direction  on  both  wings,  with  the  pitching  moment 
contributions  being  of  opposite  sign  on  the  two  wings. 


LATERAL-DIRECTIONAL 

Comparisons  between  linear  theory  calculation  and  measurement  are  shown 
In  figure  47  for  a  forward  wing  Incidence,  of  -0.27  degree  (first  test, 

{  r 

TWT-438)  and  on  figure  48  for  wc  =  -1.27  degrees  (second  test,  TWT-487). 

•  c  c 

Yawing  moment  and  side  force  derivatives  Op  and  y^  show  weak  variations  with 

angle-of-attack  and  sideslip  in  the  ranges  -4  <  a  <  +  8  and  -10  <  p  <  +  10 
d.gr.«.  Rolling  moment  dorlvative  C,  stlws  a  var(at(on  wUh  ang,a  „«  attack 

characteristic  of  a  sweptback  wing,  me  measured  derivatives  unB  and  ’y„  were 
within  ±10  percent  of  the  linear  theory  calculation  as  showrr  in  taole  V. 

The  predicted  derivative  ci|!  “  “  =  0  de9rees  »•  »»t1sfactory. 
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LONGITUDINAL  CONTROL 

The  lift  and  pitching  moment  effectiveness  of  20  percent  chord  trailing 
edge  flaps  showed  relatively  weak  variations  with  angle  of  attack  in  the  range 
-4  degrees  <  a  <  +8  degrees  and  are  lower  than  the  linear  theory  calculation  as 
shown  on  figures  49  anu  50,  on  tables  VI  and  VII  and  in  reference  8,  Measured 
increments  were  generally  in  the  range  of  60  to  90  percent  of  linear  theory 
predictions.  The  difference  between  calculation  and  measurement  is  attributed 
to  viscous  and  nonlinear  compressibility  effects. 


LATERAL  CONTROL 

The  rolling  moment  effectiveness  of  various  combinations  of  20  percent 
chord  trailing  edge  flaps  showed  weak  variations  with  angle  of  attack  in  the 
range  -4  <  a  <  +8  degrees  and  were  substantially  lower  than  the  linear 
theory  calculation  as  shown  on  figures  51,  on  table  VIII  and  in  reference  8. 
Measured  increments  were  41  to  62  percent  of  the  linear  theory  rolling 
moment  increments. 


DIRECT  SIDE  FORCE 

The  maximum  side  force  flap  configuration  is  with  the  trailing  edges  of 
all  three  flaps  deflected  toward  each  other  on  one  side  of  the  joined  wing 
and  away  from  each  other  on  the  other  side,  i.e. ,  with: 


+  S2 


*  ®*> 


L 

*R 


+  65 
+  \ 


L 


+  ^3^  =  +^F 

+  ^3„  =  -  8F 

+  S*L  .  _6A 

*  \  ‘  ^ 


L 

R 

L 

R 
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f  t  A  /i 

For  F,  =  An  =  +10  degrees  and  FR  =  A.  =  -10  degrees,  .comparisons 
between  linear  calculation  and  measurement  areL(M  =  0.90,  a  =  0  degrees): 


Linear  Theory 

Measurement 

acl 

0.0 

-0.026 

aCm 

0.0 

-0. 016 

ACy 

+0.0306 

+0.018 

A^n 

+0.0059 

+0.0055 

AC, 

+0.0043 

+0.00385 

As  with  other  flap  deflections,  the  measured  increments  are  less  than 
calculated  by  linear  theory.  Since  it  is  of  interest  to  consider  flap  deflec¬ 
tions  to  reduce  the  yawing  and  rolling  moments  to  zero  while  maintaing  a  side 
force  increment,  the  linear  theory  results  are  given  in  equation  form  as: 


C,  =  +0. 0053456c  +0.005796.  +0.0053456c 

L  r|_  ‘ r 

CM  =  +0.010606p  -0.0084356ft  +G.Q10606j: 

Cy  *  +0.000736p  -0.C00806a  -0.000736f 

Cn  -  +0. 0006356c  +0.000346.  -0.0006356c 

n  r  ^ 

C£  =  +0.001006F  +0.0007856a  -0.001006p 

L  L  R 


By  inspection  it  can  be  seen  that  the  flap  deflections 
6r  _  -6f  _  -0,5536.  _  +0.5536. 

fL  -  t-R  -  \  “  AR 

* 

give  zero  yawing  moment,  and: 

+6C  _  -6C  -  *0.7876.  _  +0.7876. 
hL  hR  \  "  aR 

give  zero  rolling  moment. 

Both  of  these  combinations  give  positive  cideforce.  The  flap  deflections 
on  the  aft  wing  are  greater  than  on  the  forward  wing  to  produce  direct  si  deforce 
with  either  zero  yawing  moment  or  zero  rolling  moment. 
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SUMMARY 


Chord  plane  linear  theory  comparisons  reported  here  have  successes  and 
deficiencies  that  are  similar  to  monoplane  correlation  experience.  That  is, 
slopes  are  satisfactorly  predicted  and  intercepts  and  control  effectiveness 
are  marginal. 


SURFACE  PRESSURE  DATA 


Pressure  data  from , the  first  wind  tunnel  test  with  the  forward  wing  at 
nominal  zero  incidence  (' wf  =  -0.27  degrees)  have  been  discussed  for  the  design 
condition  N  =  0.90  and  lift  coefficient  0.50  in  the  section  on  transonic 
redesign.  The  shock  on  the  upper  surface  of  the  forward  wing  was  predicted 
by  the  full  potential  analysis,  but  was  not  captured  by  the  transonic 
small  disturbance,  analysis.  With  a  one  degree  nose  down  change  in  the 
Incidence  of  the  forward  wing  to  'wF  *  -1.27  degrees,  the  full  potential  pre¬ 
diction  Indicated  the  shock  would  be  eliminated  on  the  forward  wing.  Surface 
pressures  measured  on  the  forward  wing  during  the  second  wind  tunnel  test  with 
this  one  degree  incidence  change  are  compared  on  figure  52  with  calculations. 
Figure  52a  shows  that  the  measured  inboard  shock  was  further  aft  than  calcu¬ 
lated.  Figure  52b  indicates  that  the  midspan  shock  was  weakened  but  not 
eliminated  by  the  incidence  change.  Figure  52c  shows  reasonable  agreement 
between  calculation  and  measurement  at  the  outboard  station  with  a  suggestion 
of  trailing  edge  separation.  The  measured  lower  surface  pressures  were 
slightly  more  positive  than  calculated  at  all  three  spanwise  stations. 
Oil  flow  photographs,  figure  53,  corroborated  the  upper  surface  shock  positions 
at  the  Inboard  and  midspan  stations. 

Surface  pressures  measured  on  the  aft  wing  are  compared  on  figure  54  with 
pressures  calculated  using  transonic  small  disturbance  analysis  with  effective 
twists  considering  downwash  from  the  forward  wing.  The  measured  pressure  data 
at  the  Inboard  station,  ligure  54a,  from  the  second  wind  tunnel  test  are 
suspect  for  reasons  which  will  be  subsequently  stated.  Also  shown  at  the 
inboard  station  are  pressure  data  from  the  first  wind  tunnel  test  at  the  same 
total  lift  coefficient.  To  produce  a  lift  coefficient  of  0.5,  the  angle-of- 
attack  was  6.0  degrees  during  the  first  wind  tunnel  test  and  6.5  degrees  during 
the  second.  The  measured  pitching  moment  was  more  negative  during  the  second 
w.nd  tunnel  test  so  that  at  the  same  lift  coefficient,  the  aft  wing  should 
carry  more  lift  and  the  forward  wing  less  lift.  The  measured  pressures  at  the 
Inboard  station  on  the  aft  wing  show  substantially  less  lift  load  being  carried 
during  the  second  test  and,  for  this  reason,  are  questionable. 

Measured  pressure  data  from  the  second  wind  tunnel  test  are  compared  with 
TSD  calculation  at  the  midspan  and  outboard  stations  on  figures  54b  and  54c  at 
an  angle-of-attack  of  6.5  degrees.  The  mid-span  of  the  aft  wing  is  carrying 
less  lift  than  calculated,  while  the  outboard  station  is  carrying  about  the 
same  lift  as  predicted.  Measurement  and  calculation  are  in  reasonable  agree¬ 
ment  for  the  more  lightly  loaded  slightly  supercritical  aft  wing.  Both 
pressure  measurements  and  oil  flow  observations  indicated  an  absence  of  shocks 
on  the  aft  wing  at  the  design  condition. 

Complete  sets  of  tabulated  and  plotted  pressure  data  are  available  in 
reference  8. 
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Section  VII 
CONCLUSIONS 


1.  The  joined  wing  supercritical  design  achieved  88  percent  of  the  fully 
turbulent  theoretical  upper  bound  L/D  of  12.97  which  neglects  compo¬ 
nent  viscous  interference  and  separation. 

2.  Weak  drag  divergence  was  achieved  at  the  design  condition.  A  swept 
shock  exists  on  the  forward  wing  upper  surface  as  a  result  of  not 
achieving  the  anticipated  twist  change  due  to  model  elasticity.  The 
impact  of  this  effect  is  estimated  to  be  0.5  in  L/D  and  was  deduced 
from  the  measurements  between  M  =  0.85  and  0.90.  Elimination  of  this 
decrement  results  in  achievement  of  92  percent  of  the  ultimate  lift/drag 
ratio. 

3.  Theoretical  turbulent  upper  bound  aerodynamic  efficiency  for  the  reference 
monoplane  and  joined  wing  are  essentially  the  same  for  the  same  mean 
sweep.  The  vortex  drag  advantage  of  the  nonplanar  joined  wing  of  the  same 
span  is  offset  by  viscous  penalties  resulting  from  higher  surface  area  for 
the  same  plan  area  and  lower  chord  Reynolds  number. 

r 

4.  Joined  wing  longitudinal  trim  is  impacted  by  elastic  effects  on  M 
resulting  from  differential  twist  between  the  front  and  aft  surfaces? 
Out  of  balance  pitching  moment  will  be  typically  reduced  at  high  dynamic 
pressure  for  unstable  balances. 

5.  Six  component  force  measurements  for  the  joined  wing  are  well  behaved  and 
near  linear  over  an  angle-of-attack  range  of  -4<  «  <  10  degrees.  Stable 
stall  is  naturally  achieved  for  the  arrangement.  Weak  changes  in  longi¬ 
tudinal  and  lateral -directional  stability  and  control  with  pitch  angle 
exist. 

6.  Chord  plane  linear  theory  comparisons  have  successes  and  deficiences 
that  are  similar  to  monoplane  correlation  experience.  That  is,  slopes 
are  satisfactorily  predicted  and  intercepts  and  control  surface 
effectiveness  are  marginal. 

7.  Approximate  and  simplified  theoretical  structural  analyses  underestimated 
the  stiffness  of  the  joined  wing  wind  tunnel  model. 
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TABLE  III  JOINED  WING  AIRFOIL 
DESIGN  CONDITIONS 

M,  -  0.90 

CL  -  0.S0 

dQn  .  0.20 

dCL 

t  -  o.os 

c 


FORWARD  PANEL 

AFT  PANEL 

ca  ■ 

0.60 

c* 

m 

0.33 

tc/2  . 

41° 

tc/z 

SB 

-48° 

Mx  - 

0.68 

MX 

m 

0.602 

tin  - 

1. 

Cij. 

m 

0.74 

-1 

0.66 

c*x 

cU 

■ 

0.747 
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Table  IV 


LONGITUDINAL  CHARACTERISTICS 

a  =  0°,  P  = 

0° 

1w 

WF 

M 

CL 

Linear  a 

Theory  Meas. 

CL 

Linear1*  0 
Theory  Meas. 

dCuJ 

Linear 

Theory 

dcL 

Meas. 

C|y| 

Linear1*  0 
Theory 

Meas. 

-0.27 

0.90 

0.0687  0.071 

0.1200  0.0736 

0.200 

0.243 

+0.0030 

-0.0041 

-1.27 

0.90 

0.0687  0.075 

0.0991  0.0445 

0.200 

0.236 

-0.0635 

-0.0569 

-1.27 

0.40 

0.0587  0.061 

0.1080  0.0563 

0.140 

0.164 

-0.0510 

-0.0439 

*As-bu11t  camber  and  twist 

Table  V 

LATERAL  DIRECTIONAL  STABILITY 
a  =  0°,  -5°  <  p  <  +5° 


r 

T  d 

LI  near  H 

Linear 

Linear  ^ 

M 

Theory 

Meas. 

Theory  v 

Meas. 

Theory 

Meas. 

-0.27 

0.90 

0.01195 

0.012 

0.0032 

0.0028 

-0.00053 

-0.00092 

-1.27 

0.90 

0.01195 

0.013 

0.0032 

0.0028 

-0.00053 

-0.00055 

-1  97 

-0.00053 

-0.00052 
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Table  VI 


LIFTING  FLAP  EFFECTIVENESS 
M  =  0.90,  a  =  0  degrees 


Flap 

Combination 

Flap 

Deflection* 

(degrees) 

1wF 

Linear 

Theory 

C|_ 

Measurement 

4+5 

+10 

-1.27 

0. 00760 

•  0.0047 

4+5 

+10 

-0.27 

0. 00760 

0.0058 

4+5 

-10 

-0.27 

0.00760 

0.0044 

1+2+4+5 

+10 

-0.27 

0.01721 

0.0144 

1+2 

+10 

-0.27 

0.00961 

0.0087 

Table  VII 

PITCHING  MOMENT  EFFECTIVENESS 
M  a  0.90,  a  «  0  degrees 


Flap 

Combination 

Flap 

Deflection* 

(degrees) 

1w 

WF 

Linear 

Theory 

^m 

Measurement 

4+5 

+10 

-1.27 

0.01264 

0.00780 

4+5 

+10 

-0.27 

0.01264 

0.00001 

4+5 

-10 

-0.27 

0.01264 

0.01070 

1+2+4+5 

+10 

-0.27 

0.00565 

0.00220 

1+2 

+10 

-0.27 

0.01829 

0.01290 

^Perpendicular  to  hinge  line 
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Table  V III 

ROLLING  MOMENT  EFFECTIVENESS 
M  =  0.90,  a  -  0  degrees 


Flap* 

Combination 

Flap 

Deflection** 

(degrees) 

V 

Linear 

Theory 

CA 

Measurement 

3R 

+10 

-0.27 

-0.000310 

-0.000127 

6R 

+10 

-0.27 

-0.000470 

-0.000218 

3R  +  6R 

+10 

-0.27 

-0.000780 

-0.000345 

3L  +  6L 

3R  +  6R 

+10 

-10 

-1.27 

+0.001560 

+0.000728 

2L  +  6L 

2R  +  6R 

+10 

-10 

-1.27 

+0.00181 

+0.001112 

*R  -  Right  side,  L  -  Left  side 

**  Perpendicular  to  hinge  line 
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Figure  1.  Reference  Monoplane  Arrangement 


Figure  2.  Joined  Wing  Arrangement 


FIGURE  4.  COMPLETED 
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Figure  5.  Effect  of  Longitudinal  Stability  on  Joined  Wing  Trimmed  Span  Load 

Efficiency  at  M  =  0.90 
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Figure  10.  Forward  Joined  Wing  Design  Pressure  Distribution  at  M  =  0.' 


Figure  12,  Yawed  Wing  Boundary  Layer  Analysis  of  Design  Pressure  Distributi 
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3)  Airfoil  at  n 
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Figure  23.  Final  Aft  Wing 
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FIGURE  24.  CCNTEIJED 


c)  Upper  Surface  Isobars 
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c)  Airfoil  With/Without  Boundary  Layer  Displacement  Thickness;  n 
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d)  Airfoil  With/Without  Boundary  Layer  Displacement  Thickness  ; 

FIGURE  26.  CONCLUDED 
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Figure  31.  Joined  Wing  Theoretical  Elastic  Twist  Increments, 

CL  =  0.50,  q  =  1120  psf 
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Figure  34.  Joined  Wing  Full  Potential  Pressure  Distribution  Comparison,  Forward 

Wing,  M=  0.90 
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Figure  36.  Joined  Wing  Transonic  Small  Disturbance  Pressure  Distribution 

Comparison,  Aft  Wing,  M  *  0.90,  CL  =  0.50 
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Figure  45.  Reference  Monoplane  Design  Pressure  Distribution  at  M  ~  0.90, 

CL  -  0.50 
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Drag  Ratio  Comparison  at  K  -  0.90, 
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<: inure  52.  Comparison  of  Measured  and  Full  Potential  Forward  Wing  Pressure 

Distributions  at  M  *  0.90,  =  0.50 
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FIGURE  52  CONTINUED 
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Appendix  A 

AERODYNAMIC  METHODOLOGY 
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LINEAR  DESIGN  AND  ANALYSIS 


The  initial  design  and  analysis  tool  used  is  the  unified  distributed 
panel  method  (reference  2),  The  linearized  small  'disturbance  equation 

(1  -  M2  )  d>  +  6  H  =  0 
v  *xx  yyy  yzz 


is  numerically  solved  for  arbitrary  configurations  by  a  superposition  of 
chord  plane  and  body  singularities  in  conjunction  with  an  interference  shell 
to  account  for  surface-body  inductions.  Surface  lifting  effects  are 
represented  by  constant-pressure  panels.  Surface  thickness  effects  are 
represented  by  chordwise  linearly  varying  source  panels.  The  analysis  is 
applicable  to  multiple-surface  nonplanar  configurations  at  subsonic  or 
supersonic  speeds.  Lorgitudinal  and  lateral-directional  forces  and  rotary 
derivative  characterises  are  evaluated. 

Inverse  lifting  solutions  are  obtained  by  specifying  the  net  loads  at 
the  panel  centroids  and  determining  the  boundary  conditions  by  matrix  multi¬ 
plication.  Problems  where  boundary  conditions  are  given  in  one  region  and 
net  loads  are  prescribed  in  the  remainder  can  also  be  treated.  Wing- on- wing 
inductions  in  the  presence  of  the  body  are  established  by  this  analysis  for 
subsequent  use  in  the  transonic  codes  to  approximately  account  for  surface 
interactions  (by  use  of  an  effective  twist  and  camber  philosophy)  not  directly 
modeled. 

Linear  analysis  is  also  used  in  conjunction  with  structural  influence 
coefficients  tc  estimate  static  aeroelastic  deflections  of  the  model  wing 
panels. 

In  the  design  mode,  lifting  surface  theory  is  used  to  minimize  (reference 
3)  the  fa^-field  ( i . e . Trefftz  plane)  vortex  drag  and  define  the  optimum 
surface  span  loads  and  candidate  twist  distributions. 


SUPERSONIC  PRESSURE  DRAG 

The.  tec unique  used  to  evaluate  wave  drag  is  the  supersonic  area  rule. 
Th«  spatial  singularities  which  are  a,  solution  to  the  linearized  equation  of 
motion  are  reduced  to  a  series  of  equivalent  linear  distributions  by  applica¬ 
tion  of  the  cutting  plane  concept.  The  3-D  distribution  is  surveyed  longi¬ 
tudinally  at  fixed  roll  angles  by  an  oblique  plane.  The  drag  is  then  evaluated 
with  slender  body  theory. 


A-l 


1 

4t?S 


2tt  L  L 


ffl  *■<*!■  e> 


S"(x2,0)  1 nCxj  -  x2]  dx1dx2d& 


where  S  is  the  oblique  projected  cross-sectional  area  cut  by  the  system  of 
Mach  planes. 


TRANSONIC  DESIGN  AND  ANALYSIS 


The  methods  used  to  derive  supercritical  forward  and  aft  wing  contours 
in  the  presence  of  each  other  and  the  body  are  described  in  the  following 
paragraphs.  Experience  with  transonic  methods  indicates  that  the  initial 
supercritical  contours  be  defined  with  the  transonic  small  disturbance  formula¬ 
tion  because  of  its  design  capability.  The  result  Is  subsequently  evaluated 
and  refined  using  a  full -potential  analysis. 


TRANSONIC  SMALL  DISTURBANCE  THEORY 

TSD  theories  are  extensions  of  the  classical  small  disturbance  formula¬ 
tion.  The  modifications  involve  the  retention  of  higher  order  terms,  princi¬ 
pally  to  capture  swept  shocks.  The  governing  equation  Is  solved  by  a  mixed 
elliptic- hyperbolic,  relaxation  algorithm  applied  to  a  finite-difference 
approximation.  The  boundary  conditions  are  linearized  and  applied  on  the 
surface  chord  plane  and  fuselage  shell. 

The  wing-body  code  (reference  Sr)  uses  an  embedded  fine-grid  structure 
with  sweep  and  taper  transformations  in  addition  to  a  crude  Castesian  outer 
grid.  This  feature  increases  computational  efficiency  by  a  factor  of  five 
over  earlier  small  disturbance  solutions  of  the  same  resolution.  The 
equation  used  for  computing  transonic  flow  is: 
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The  underlined  terms  are  those  added  to  the  classical  small  disturbance 
formulation. 

s 

Rockwell's  wing-body  inverse  code  (reference  6)  is  used  to  define  the  wing 
contours  for  specified  design  pressure  distributions. 


FULL  POTENTIAL  ANALYSIS 

The  Jameson  FL028  transonic  wing-body  (reference  9)  code  solves  the  exact 
potential  equation. 


^x.  K  v  ~  *  0 

X1  xj  xixj  XjX. 

U  j  -  1,2P3 


using  a  conservative  rotated  differencing  scheme  for  a  surface  fitted  grid 
system  obtained  from  a  series  of  analytic  mappings. 


WING  BOUNDARY- LAYER  ANALYSIS 

A  yawad  wing  integral  laminar/turbulent  boundary  layer  analysis  (reference 
4)  is  used  to  assess  candidate  sectional  design  pressure  distributions  from 
an  attached  flow  viewpoint  at  the  test  condition.  The  solution  is  also  used  to 
estimate  transition  point  location  considering  cross  flow  vorticity  due  to 
sweep  and  streamline  Instability.  Finally,  the  evaluation  Is  used  to  estimate 
the  magnitude  of  the  displacement  thickness  correction  that  will  be  required 
to  eliminate  weak  viscous  boundary  layer  Interaction  effects  by  undercutting 
the  design  section  geometry. 

The  three  dimensional  laminar  and  turbulent  boundary  layer  analysis 
developed  by  Nash  (reference  7)  is  applicable  to  finite  swept  wings.  An 

orthogonal  surface  coordinate  system  is  used,  and  the  equations  are  solved 

with  an  efficient  difference  scheme.  The  analysis  is  used  to  verify  that  the 
final  win<j  design  pressure  distributions  produce  attached  flow  at  the  test 
Reynolds  number  and  estimate  potential  separation  problems  due  to  adverse 

pressure  gradients  and  shock  boundary- layer  interactions  which  can  not  be 

fully  controlled.  For  the  design  case  the  inviscid/vicous  interaction  is 
minimized  by  removing  the  three  dimensional  displacement  thickness  from  the 
wing  sections. 
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Appendix  B 

VISCOUS  DRAG  ANALYSIS 


Standard  techniques  (reference  2)  for  the  evaluation  of  turbulent  flat 
plate  skin  friction  at  compressible  speeds  are  used  4in  conjunction  with  a 
component  buildup  approach.  Fuselage,  empennage,  and  wing  fairing  component 
thickness  effects  were  approximated  using  experimental  data  correlations.  An 
Integral  infinite  yawed  boundary  layer  analysis  (reference  4)  is  used  to 
correct  the  flat  plate  wing  analysis  for  the  effect  of  thickness,  lift  and 
sweep  as  follows 


Z  "f.  £l 
.  ’  Sref 


»  -c  L 


The  mean  chord  sectional  profile  drag  LD  is  evaluated  by  relating  the 
wake  momentum  thickness  at  Infinity  to  the  momentum  thickness  at  the  trailing 
edge  using  a  Squire-Young  analysis  appropriately  extended  for  the  effect  of 
sweep.  That  is 


'dp  *  E 


where 


per 


<H  +  H  +  4) 


--  T 


|"HTE  + 

L  - - 


U'TE  “  Vt  [\A>2  -  slAc/Z] 


TE  u  TE'  Moo 

■  uw  " 
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and  Hto  »  1  +  (y-l )  m2 
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JOINED  WING  SECTIONAL  BOUNDARY  LAYER  AND 
PRGFILE  DRAG  CHARACTERISTICS 
M  =  0.90,  RN/Vx  =  7  x  106 
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NOMENCLATURE 

r 

dp  Sectional  flat  plate  skin  friction  drag  coefficient 
r 

dp  Sectional  profile  drag  coefficient 

r 

Dp  Strip  flat  plate  skin  friction  drag  coefficient 
r 

Dp  Profile  drag  coefficient 

Cp  Average  flat  plate  skin  friction  coefficient 

c  Mean  aerodynamic  chord 

a 

H  Boundary  layer  form  factor.  6 

M  Mach  number 

d  Total  velocity 

Sref  Reference  area 

4$  Stri p  area 

T  Static  temperature 

y  Ratio  of  specific  heat,  1.4 

6*  Boundary  layer  displacement  thickness  in  potential  streamline 

direction 

0,,  Boundary  layer  momentum  thickness  in  potential  streamline 
direction 

A  Wing  sweep 


SUBSCRIPTS 


%.  Half  chord 

e  Edge  of  boundary  layer 

1  Wing  strip  1 

TE  Trailing  edge 

»  Free  stream 
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Appendix  C 

MODEL  DESIGN  AND  CONSTRUCTION 


STRESS  AND  FLUTTER  ANALYSES 

Air  loads  for  the  l/10th  scale  joined  wing  arrangement  of  figure  2 
exceeded  available  force  balance  capacity  so  the  model  scale  was  reduced  to 
l/15th.  This  changed  the  overall  span  of  both  the  reference  monoplane  and  the 
joined  wing  from  54  Inches  to  36  Inches.  Corresponding  reductions  in 
reference  area  were  from  911.25  square  Inches  down  to  405.00  square  Inches  and 
the  reference  mean  aerodynamic  chord  (m.a.c)  from  18.51  inches  to  12.34 
inches.  A  limit  condition  of  Mach  number  0.9,  angle  of  attack  10  degrees, 
lift  coefficient  0.785  and  dynamic  pressure  1120  pounds/square  foot  gave  a 
limit  11ft  lead  on  the  model  of  2500  pounds. 

The  redundant  structure  of  the  joined  wing  model  required  special  stress, 
flutter  and  deflection  analyses.  Simple  beam  model  stress  analyses  indicated 
excessive  stress  levels  near  the  wing  tips  where  the  two  wings  joined.  To 
relieve  this  stress  problem  and  to  provide  adequate  thickness  for  installa¬ 
tion  of  pressure  orifices  near  the  tips,  the  thickness  ratio  at  the  wing  tips 
was  increased  from  0.05  up  to  0.07,  varying  linearly  to  0.05  at  the  model 
midspan. 

A  finite  beam  element  flutter  theoretical  analysis  was  conducted, 
reference  5,  and  It  was  concluded  that  wing  tip  fairings  of  either  aluminum 
or  steel,  must  have  a  width  of  at  least  5  percent  the  joined  wing  semispan, 
to  provide  good  bearing  surfaces  and  rigid  wing-to- fairing  attachment  for 
adequate  flutter  margin. 


AERGELASTIC  ANALYSIS 

The  structural  finite  beam  element  model  created  for  the  flutter  analysis 
was  adapted  to  the  calculation  of  static  aeroelastic  deflections  by  eliminating 
the  sting  and  balance  flexibilities,  first  two.  The  spanwise  design  load 
distributions  on  the  two  wings  were  segmented  to  give  a  series  of  concentrated 
loads  at  the  control  points  of .  the  structural  finite  element  model.  The 
resulting  angle-of-nttack  change  due  to  flexibility  (Aoflex)  is  presented  on 
figure  Cl  and  was  subtracted  from  the  design  twist  distribution  to  give  the 
jig  shape  for  model  construction.  The  intent  was  at  the  design  point,  the 
model  would  flex  Into  the  designed  twist  distribution. 


MODEL  LOFT  GEOMETRY 

Lofted  model  drawings  are  presented  In  flgurer  C2  through  C8.  Figure  C2 
presents  the  ronMgurtion  two  view  and  figure  C3  the  forward  wing  detail, 
Figures  C4  and  C5  show  the  aft  wing  detail  and  rigging  respectively.  Figure  Ce 
presents  the  vertical  tail  and  wing  tip  fairing  details. 
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MODEL  FABRICATION 

Templates  were  constructed  from  the  lofted  model  drawings.  The  model 
was  machined  to  template  contours.  Both  wings  and  the  aft-wing-to-verticaL 
attach  fitting  were  built  of  high  strength  steel.  The  fuselage  and  wing  tip 
center  sections  were  built  of  regular  strength  steel.  The  nose  fairing  of  the 
fuselage  and  the  nose  and  tall  portions  at  the  wing  tip  center  fairings  were 
built  of  aluminum.  During  model  assembly,  it  became  evident  the  original 
wing  tip  fairings  of  circular  cross  section  lacked  sufficient  depth  for  an 
adequately  stiff  wing  tip  joining.  The  circular  cross  section  'of 
0.90  Inch  diameter  was  changed  to  a  rounded  rectangular  cross  section  of 
0.95  inch  width  and  1.21  inch  depth.  The  forebody  lengths  of  the  wing  tip 
fairings  were  Increased  from  1.79  Inches  up  to  2.40  Inches  and  the  afterbody 
lengths  from  2.69  Inches  up  to  3.60  Inches  In  order  to  maintain  an  equivalent 
forebody  fineness  ratio  of  2.0  and  afterbody  fineness  ratio  of  3.0. 


CONTOUR  CHECKS 

Static  pressure  orifice's  were  installed  in  the  right  forward  and  aft 
joined  wings  at  spanwlse  at  stations  2Y/b  of  0.18,  0.51  and  0.87.  After 
model  assembly,  the  twist  angles  of  the  wings  were  measured  at  the  spanwlse 
stations  of  pressure  measurement  and  found  to  be  within  0.05  degrees  of  the 
lofted  twljt  on  the  forward  wing  and  within  0.24  degrees  of  the  lofted  twist 
on  the  aft  wing.  The  upper  surface  contour  of  the  forward  wing  was  measured 
with  a  cordax  machine  and  found  to  agree  with  the  template  contour  within 
loss  than  a  line  width  on  the  lofted  drawing. 


AEROELASTiC  MEASUREMENT 

Strain  gauges  to  measure  bending  moment  and  torque  were  Installed  on  the 
lower  surface  of  the  left  forward  wing  at  spanwlse  stations  8  and  12  inches 
off  the  model  centerline.  These  gauges  wore  calibrated  by  applying  known 
loads  to  the  Isolated  forward  wing  mounted  in  a  test  rig.  Then  with  the  model 
completely  assembled,  known  100  pounds  down  loads  were  applied  11.08  inches  off 
the  model  centerline  at  x/c  ®  0.45  perptndicula;  to  the  planas  of  the  wings  as 
shown  on  figure  C7.  Deflections  and  twists  of  the  wings  relative  to  the 
fuselage  were  measured  at  the  midspan  and  it  the  tip.  A  comparison  of  measured 
and  calculated  deflections  and  twists  Is  presented  on  table  Cl.  The  "as  built" 
model  was  substantially  stiffer  than  calculated  by  the  finite-beam-element 
theoretical  calculations  and  the  bending  moments  and  torques  measured  during 
the  known  loading  and  during  wind  tunnel  test  were  substantially  lower  than 
calculated  (figures  C8  and  C9).  These  differences  were  attributed  to 
Inaccuracies  in  the  calculated  stiffness  coefficients  and  to  Inadequacies  In 
the  theoretical  formulation  (for  example,  longitudinal  deflections  were 
neglected).  Substantial  hysteresis  was  observed  in  the  measured  bending 
moments  and  torques  due,  apparently,  to  adjustments  in  the  tightness  of  tho 
screw  and  bolt  fasteners  at  the  model  joints. 
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FIGURE  C7  WING  DEFLECTION  CALIBRATION 


BENDING  MOMENT  DISTRLft  PON  AT  LIMIT  LOADING  CONDITION 


